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ABSTRACT 

Recent investigations of the WD + MS channel of Type la supernovae (SNe la) imply 
that this channel may be the main contribution to the old population (>1 Gyr) of SNe 
la. In the WD + MS channel, the WD could accrete material from a main-sequence or 
a slightly evolved star until it reaches the Chandrasekhar mass limit. The companions 
in this channel would survive after SN explosion and show distinguishing properties. 
In this Letter, based on SN la production regions of the WD + MS channel and three 
formation channels of WD + MS systems, we performed a detailed binary population 
synthesis study to obtain the properties of the surviving companions. The properties 
can be verified by future observations. We find that the surviving companions of the 
old SNe la have a low mass, which provides a possible way to explain the formation 
of the population of single low-mass WDs (<0.45M Q ). 
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1 INTRODUCTION 

Type la supernovae (SNe la) appear to be good cosmologi- 
cal distance indicators owing to their high luminosities and 
remarkable uniformity, and have been applied successfully in 
determining cosmological parameters (e.g. fi and A; Riess et 
al. 1998; Perlmutter et al. 1999). However, several key issues 
related to the nature of their progenitors and the physics 
of the explosion mechanisms are still not well understood 
(Hillebrandt & Niemeyer 2000; Wang et al. 2008; Podsiad- 
lowski 2010), and no SN la progenitor system before the 
explosion has been conclusively identified. These uncertain- 
ties may raise doubts about the distance calibration which 
is purely empirical and based on the SN la sample of the 
low red-shift Universe (z < 0.05; Phillips 1993). 

It is widely accepted that SNe la are thermonuclear 
explosions of carbon-oxygen white dwarfs (CO WDs) in 
binaries (for the review see Nomoto, Iwamoto & Kishi- 
moto 1997). Over the past few decades, two families of SN 
la progenitor models have been proposed, i.e. the double- 
degenerate (DD) and single-degenerate (SD) models. In the 
DD model, two CO WDs with a total mass larger than the 
Chandrasekhar (Ch) mass limit may coalesce and then ex- 
plode as an SN la (Iben & Tutukov 1984; Webbink 1984; Han 
1998). Although it is suggested that the DD model likely 
leads to an accretion-induced collapse rather than to an SN 
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la (Nomoto & Iben 1985) , it is still too early to exclude the 
model as it may contribute to a few SNe la (Howell et al. 
2006; Pakmor et al. 2010). In the SD model, the compan- 
ion is probably a main-sequence (MS) or a slightly evolved 
star (WD + MS channel), or a red-giant (RG) star (WD 
+ RG channel) (e.g. Whelan & Iben 1973; Hachisu, Kato & 
Nomoto 1996; Li & van den Heuvel 1997; Yungelson & Livio 
1998; Langer et al. 2000; Han & Podsiadlowski 2004, 2006; 
Chen & Li 2007; Lii et al. 2009; Wang, Li & Han 2010, here- 
after WLH10; Meng & Yang 2010a; Wang & Han 2010a). 
Meanwhile, a CO WD may also accrete material from a He 
star to increase its mass to the Ch mass limit (WD + He star 
channel; Solheim & Yungelson 2005; Wang et al. 2009a,b; 
Ruiter, Belczynski & Fryer 2009; Wang & Han 2010b). An 
explosion following the merger of two WDs would leave no 
remnant, while the companion in the SD model would sur- 
vive and potentially be identifiable. A surviving companion 
in the SD model would evolve to a WD finally, and Hansen 
(2003) suggested that the SD model could potentially ex- 
plain the properties of halo WDs, e.g. their space density 
and ages. Note that there has been no conclusive proof yet 
that any individual object is the surviving companion of an 
SN la. It would be a promising method to test SN la pro- 
genitor models by identifying their surviving companions. 

By considering the effect of the thermal-viscous in- 
stability of accretion disk on the evolution of WD bina- 
ries, WLH10 recently enlarged the regions of the WD + 
MS channel for producing SNe la. According to a de- 
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tailed binary population synthesis (BPS) approach, WLH10 
found that this channel is effective for producing SNe la 
(~1.8 x 10~ 3 yr _1 in the Galaxy). This study also implied 
that the WD + MS channel may be the main contribution 
to the old population (>1 Gyr) of SNe la. The companions 
in this channel would survive and show distinguishing prop- 
erties. 

At present, the existence of a population of single low- 
mass (<O.45M ) WDs (LMWDs) is supported by some 
observations (e.g. Kilic et al. 2007a). (1) The first two 
LMWDs have been implied by the work of Marsh et al. 
(1995). They carried out radial velocity measurements of 7 
LMWDs, and found that two of the WDs, WD1353+409 
and WD1614+136 do not show any radial velocity varia- 
tions. (2) Maxted et al. (2000) also presented 14 LMWDs 
with no detectable radial velocity variations from a larger 
radial velocity survey of 71 WDs. (3) A recent analysis of 
348 H atmosphere WDs from the Palomar Green (PG) sur- 
vey has revealed 30 LMWDs (Liebert et al. 2005). Kilic et al. 
(2007a) concluded that 47% of the PG LMWDs searched for 
companions seem to be single. (4) The ESO SN la Progen- 
itor Survey (SPY) searched for radial velocity variations in 
more than a thousand WDs using the Very Large Telescope 
(Napiwotzki et al. 2001). Kilic et al. (2007a) also concluded 
that 15 of 26 LMWDs discovered in the SPY project do not 
show any radial velocity variations, corresponding to a single 
LMWD fraction of 58% (also Napiwotzki et al. 2007). 

The formation of the single LMWDs is still unclear. 
It is suggested that the single LMWDs could be produced 
by single old metal-rich stars which experience significant 
mass-loss prior to the He flash (Kalirai et al. 2007; Kilic 
et al. 2007a). However, within the age of the Universe, it 
is almost certainly impossible for the single stars to pro- 
duce WDs with mass close to 0.3 M Q , or even some ex- 
tremely LMWDs with mass as low as 0.2 Mq (e.g. Kilic et 
al. 2007a,b; Justham et al. 2009). Furthermore, the study of 
initial-final mass relation for stars by Han et al. (1994) im- 
plied that only LMWDs with masses >0.4Mq might be pro- 
duced from such a single-star channel, even at high metal- 
licity (Meng, Chen & Han 2008). Thus, it would be difficult 
to conclude that single stars can produce LMWDs. Justham 
et al. (2009) recently inferred an attractive formation chan- 
nel for the single LMWDs, which could have been formed 
in binaries where their companions have exploded as SNe 
la. Note that Nelemans & Tauris (1998) also proposed an 
alternative scenario to form single LMWDs from a solar-like 
star accompanied by a massive planet, or a brown dwarf, in 
a relatively close orbit. 

Han (2008) obtained many properties of the surviving 
companions of the SNe la with intermediate delay times 
(100 Myr— 1 Gyr; the delay times of SNe la are defined as 
the time intervals between the star formation and SN la 
explosion). Wang & Han (2009) studied the properties of the 
companions of the SNe la with short delay times (<100 Myr) 
from the WD + He star channel. The properties can be 
verified by future observations (e.g. the masses, the spatial 
velocities, the effective temperatures, the luminosities, the 
surface gravities, etc; referring to Wang & Han 2009). The 
purpose of this Letter is to investigate the properties of the 
surviving companions of the old SNe la from the WD + 
MS channel inferred in WLH10, and to explore whether the 
surviving companions, later in their evolution, could explain 



the existing population of single LMWDs. In Section 2, we 
describe the BPS method for obtaining the properties of the 
companions. The BPS results and discussion are given in 
Section 3. 



2 BINARY POPULATION SYNTHESIS 

In the WD + MS channel, the progenitor of an SN la is either 
a close WD + MS or WD + subgiant system, which has most 
likely emerged from the common envelope (CE) evolution of 
a binary involving a giant star. The CE ejection is still an 
open problem. Similar to the work of Wang et al. (2009b), we 
also use the standard energy equations (Webbink 1984) to 
calculate the output of the CE phase. For this prescription of 
the CE ejection, there are two highly uncertain parameters, 
i.e. a cc and A, where a cc is the CE ejection efficiency, and A 
is a structure parameter that depends on the evolutionary 
stage of the donor. As in previous studies, we combine a cc 
and A into one free parameter q cc A, and set it to be 0.5 (for 
details see Wang et al. 2009b). 

To obtain the distributions of properties of the surviv- 
ing companions, we performed a Monte Carlo simulation 
in the BPS study. In the simulation, by using the Hurley's 
rapid binary evolution code (Hurley et al. 2000, 2002), we 
followed the evolution of 1 x 10 7 sample binaries from the 
star formation to the formation of the WD + MS systems 
according to three evolutionary channels (i.e. the He star 
channel, the EAGB channel and the TPAGB channel; the 
three channels here all refer to the WD progenitor's evo- 
lutionary phase when it encounters a Roche lobe overflow 
(RLOF) event, for details see WLH10). If a binary evolves 
to a WD + MS system, and if the system, at the onset of the 
RLOF phase, is located in the SN la production regions (see 
Fig. 6a of WLH10) in the plane of (log P 1 , Mj) for its M^ D , 
where P 1 , M\ and M^d are, respectively, the orbital period, 
the secondary's mass, and the WD's mass of the WD + MS 
system at the onset of the RLOF, we assume that an SN la 
is resulted, and the properties of the WD + MS system at 
the moment of SN explosion are obtained by interpolation in 
the three-dimensional grid (Mjy D , logP 1 ) of the ~2400 
close WD + MS systems calculated in WLH10. 

In the BPS study, the primordial binary samples are 
generated in the Monte Carlo way. We adopted the following 
input for the simulation (e.g. Wang et al. 2009b, 2010). 

(1) The initial mass function (IMF) of Miller & Scalo 
(1979) is adopted. 

(2) We take a constant mass-ratio (q') distribution 
(Goldberg & Mazeh 1994), 

n(q') = 1, < q' ^ 1, (1) 

where q' = M%/Mf. 

(3) We assume that all stars are members of binaries 
and that the distribution of separations is constant in log a 
for wide binaries, where a is separation and falls off smoothly 
at a small separation 

a-n(a) = ( a ^ a / a °^> a < a °< ( 2 ) 
[ asep, ao < a < Ol, 

where a scp « 0.07, a = lOfi ,oi = 5.75 xlO 6 R Q = 0.13 pc 
and m w 1.2. This distribution implies that the numbers of 
wide binaries per logarithmic interval are equal, and that 
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Figure 1. Distribution of properties of the companions in the 
plane of (V OI b, Mj N ) at the current epoch, where V OI h is the or- 
bital velocity and M| N the mass at the moment of SN explosion. 




about 50% of stellar systems have orbital periods less than 
100 yr (Han, Podsiadlowski & Eggleton 1995). 

(4) A circular orbit is assumed for all binaries. The or- 
bits of semidetached binaries are generally circularized by 
tidal forces on a timescale which is much smaller than the 
nuclear timescale. 

(5) The star-formation rate (SFR) is taken to be con- 
stant over the past 15 Gyr. 



3 RESULTS AND DISCUSSION 

The simulation gives current-epoch distributions of many 
properties of companions at the moment of SN explosion, 
e.g. the masses, the orbital periods, the orbital separations, 
the orbital velocities, the effective temperatures, the lumi- 
nosities, the surface gravities, the surface abundances, the 
mass-transfer rates, the mass-loss rates of the optically thick 
stellar winds, etc. The simulation also shows the initial pa- 
rameters of the primordial binaries and the WD + MS sys- 
tems that lead to SNe la. Figs 1-5 are selected distributions 
that may be helpful for identifying the surviving compan- 
ions. 

Fig. 1 shows the distributions of the masses and the 
orbital velocities of companions at the moment of SN explo- 
sion. In the figure, the companion has an orbital velocity of 
~80— 220km/s for a corresponding mass of ~0.2— 1.8 Mq at 
the moment of SN explosion. The distributions correspond 
to the moment of SN explosion. The physical quantities de- 
picted in the figure were calculated based on the assump- 
tion that the companion (e.g. at the time when the WD 
explodes) has not yet been affected by the explosion itself. 
Since SN explosion is expected to have a direct impact upon 
the surviving companion, one would expect that the distri- 
butions would look somewhat modified at a later time. For 
example, the SN ejecta will interact with its companion. 
The companions will be stripped of some mass (see the next 
paragraph) and receive a kick velocity that is perpendicular 
to the orbital velocity. Marietta, Burrows & Fryxell (2000) 
presented several high-resolution two-dimensional numerical 
simulations of the impact of SN la explosion with compan- 
ions. The study implied that this impact makes the compan- 
ion in the WD + MS channel receive a kick of 49— 86km/s. 



Figure 2. Bimodal distribution of the mass for the companions 
at the moment of SN explosion. 



With detailed stellar models and realistic separations that 
were obtained from binary evolution, Meng, Chen & Han 
(2007) obtained a similar kick velocity 30— 90km/s. Thus, 
a surviving companion has a space velocity larger by ~10% 
than that in Fig. 1. 

We find that the distribution of the mass for the com- 
panions is bimodal (see Fig. 2). The left peak results from 
the companions of SNe la with long delay times, while the 
right from the companions of SNe la with intermediate delay 
times. For SNe la with intermediate delay times, SN la ex- 
plosions occur when the companion is a MS or a slightly 
evolved star. The study of Marietta, Burrows & Fryxell 
(2000) implied that the impact of the SN explosion makes 
the MS or the slightly evolved companion lose a mass of 
0.15-0.17 M Q , while Meng, Chen & Han (2007) obtained 
a lower 'stripped mass' ~0.03— 0.13M Q . A surviving com- 
panion of SNe la with intermediate delay times therefore 
has a mass lower by ~0.1 Mq. However, for the SNe la with 
long delay times, SN la explosions occur when the compan- 
ion evolves to the RG stage. Marietta, Burrows & Fryxell 
(2000) found that an RG donor will lose almost its entire 
envelope (96%— 98%) owing to the impact of the SN explo- 
sion and leave only the core of the star. Thus, the surviving 
companions of SNe la with long delay times will contribute 
to the population of single LMWDs. 

WLH10 inferred that the WD + RG channel of SNe la, 
in which the WD begins to accrete material from an RG 
star, has a long delay time from the star formation to SN 
explosion due to the low initial mass (<1.5Mq) of the RG 
donors. So, all SNe la from this channel also contribute to 
the old populations of SNe la, and the companions have 
a low mass at the moment of SN explosion (see Fig. 2 of 
Wang & Han 2010a), which also provides a possible path- 
way for the formation of the single LMWDs. However, the 
Galactic SN la birthrate from the WD + RG channel is low 
(~3x 10 yr - ; WLH10) compared with observations, so in 
this study we mainly focus on the investigation of surviving 
companions from the WD + MS channel, which may be the 
main contribution to the population of single LMWDs (note 
that the theoretical birthrate from the WD + RG channel 
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Figure 3. Same as Fig. 1, but in the plane of (logT ett -, logp), 
where T e g is the effective temperature of the companions at the 
moment of SN explosion and log g the surface gravity. The error 
bars denote the location of Tycho G (Ruiz-Lapucnte et al. 2004). 
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Figure 4. Same as Fig. 1, but in the plane of (logP SN , M| N ), 
where P SN is the orbital period at the moment of SN explosion. 
The filled asterisk indicates the position of a recurrent nova (U 
Sco). 



still contains many potential uncertainties). Here, we em- 
phasize that the surviving companions of the old SNe la 
from the WD + MS and WD + RG channels provide a pos- 
sible way for the formation of the population of the single 
LMWDs. We also suggest that the observed single LMWDs 
may provide evidence that at least some SN la explosions 
have occurred with non-degenerate donors (such as MS or 
RG donors). 

The timescale from SN explosion to the formation of a 
WD from the surviving companion is important in future ob- 
servational tests, which is related to the evolutionary phase 
of the surviving companion. If the companion is in the MS 
stage at the moment of SN explosion, the timescale mainly 
depends on the nuclear burning lifetime of the star. We did a 
test for the evolution of a typical MS surviving companion, 
which has the mass of 1.04 Mq and with central H abun- 
dance 0.26 (its remaining MS lifetime is about 840 Myr, and 
its post MS lifetime is about 1300 Myr). Thus, on average 
the nuclear burning lifetime of the MS surviving compan- 
ion is about 2 Gyr. If the companion is in the RG stage at 
the moment of SN explosion, the timescale is mainly de- 
cided by the lifetime of the H-shell burning. We also did a 
test for this case (e.g. a 0.42 Mq RG companion star with 
0.17 Mq He-core, the lifetime of the H-shell burning is about 
180 Myr. Here, we ignore the thermal equilibrium timescale 
from a He-core to a WD , since it is short ~10 6 yr compared 
with the timescale of the H-shell burning). Therefore, a low 
mass companion in the RG stage will evolve to a WD more 
quickly than the massive companion in the MS stage. 

Fig. 3 represents the distributions of the effective tem- 
peratures and the surface gravities of the companions at the 
moment of SN explosion. Tycho G was taken as the surviv- 
ing companion of Tycho's SN by Ruiz-Lapuente et al. (2004). 
It has a space velocity of 136km/s, more than 3 times the 
mean velocity of the stars in the vicinity. Its surface gravity 
is log (g/cms~ 2 ) = 3.5 ±0.5, while the effective temperature 
is T e ff = 5750 ± 250K. The parameters are compatible with 
Figs 1 and 3. However, Fuhrmann (2005) argued that Tycho 
G might be a Milky way thick-disk star that is coincidentally 
passing the vicinity of the remnant of Tycho's SN. Ihara et 
al. (2007) recently also argued that Tycho G may not be the 
companion of Tycho's SN, as the star did not show any spe- 



cial properties in its spectrum. The surviving companions of 
SNe la would expect to be contaminated by SN ejecta and 
show some special characteristics (e.g. Marietta, Burrows 
& Fryxell 2000). Thus, whether Tycho G is the surviving 
companion of Tycho's SN is still quite debatable. 

Fig. 4 shows the distributions of orbital periods and 
secondary masses of the WD + MS systems at the moment 
of SN explosion. The orbital periods and secondary masses 
at this moment are basic input parameters when one simu- 
lates the interaction between SN ejecta and its companion. 
This figure may also help us to verify whether some WD + 
MS systems observed could explode as SNe la. A recurrent 
nova (U Sco) is indicated by a filled asterisk in the figure, 
in which the WD mass is about 1.37 Mq and its compan- 
ion is a 1.5 M MS star (Hachisu et al. 2000). The orbital 
period of the binary is 1.23 d (Schaefer & Ringwald 1995). 
Hachisu et al. (2000) concluded that the WD could increase 
its mass until an SN la explosion. If the WD explodes as 
an SN la eventually, the companion would have a slightly 
smaller mass. The binary orbital period will first decrease 
and then increase if the mass-ratio reverses. Then, its final 
position in the figure will move to a lower mass than its 
present one, and may enter into the most probable area for 
producing SNe la. Thus, U Sco is likely to explode as an SN 
la (see also Meng & Yang 2010b). 

If we assume that the companions co-rotate with their 
orbits, we can obtain the distributions of their equatorial 
rotational velocities (see Fig. 5). We see that the surviving 
companions are fast rotators, so their spectral lines should 
be broadened noticeably. The rotational velocity of com- 
panions from the WD + MS channel is in the range of 
~10— 140km/s, which is lower than that from the WD + 
He star channel (~120-380 km/s; Wang & Han 2009). This 
is because the WD + He star channel has shorter orbit pe- 
riods than that of the WD + MS channel at the moment of 
SN explosion. 

The simulation in this Letter was made with a co A = 0.5. 
If we adopt a higher value for a co A (e.g. 1.5), the birthrate of 
SNe la will be lower than the case of a ce A = 0.5 (the binaries 
emerged from the CE ejections tend to have slightly closer 
orbits for c* cc A = 0.5 and are more likely to be located in 
the SN la production region). In addition, a high value of 
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Figure 5. Same as Fig. 1, but in the plane of (V ro t, M§ N ), where 
Vrot is the equatorial rotational velocity of the companions at the 
moment of SN explosion. 

a cc A leads to a systematically later explosion time, i.e. the 
delay time from the star formation to SN explosion will be 
longer. This is because a high value of a cc A leads to wider 
WD + MS systems, and, as a consequence, it takes a longer 
time for the companion to evolve to fill its Roche lobe. For 
a high value of a ce X (e.g. 1.5), the minimum delay time is 
~360Myr, while the value is ~280Myr for a cc A = 0.5. 
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